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PREFACE 

A  modern  high  speed  digital  computer  uses  many  vacuum  tubes  in 
bistable  circuits.  These  circuits  and  hence  the  computer,  are  susep- 
tible  to  the  failures  of  the  vacuum  tubes.  Although  considerable 
progress  has  been  made  towards  improving  the  reliability  of  vacuum 
tubes ,  they  do  not  approach  that  of  the  later -developed  "solid- 
state"  devices.  In  addition  to  greater  reliability  offered  by  these 
devices,  the  small  size  and  low  power  consumption  is  particularly 
enticing  where  large  numbers  of  bistable  circuits  are  needed.  In 
particular,  the  development  of  the  ferrites  and  titanates  for  high 
permeability  and  dielectric  constants  introduced  non-linear  elements 
in  a  component  size  that  surpasses  the  highly  lauded  transistor.  As 
with  other  new  circuit  elements  much  is  left  to  be  done,  especially 
in  improvement  of  the  materials  involved.  However,  the  first  step 
was  taken  by  C.  F.  Spitzer  (TJ   in  the  analysis  of  the  two-branch, 
bistable,  series  resonant  circuit,  using  a  non-linear  inductance. 
The  next  move  is  the  analysis  of  the  non-linear  capacitor  circuit 
which  this  paper  purports  to  do.  The  next  obvious  step  would  be  the 
inclusion  of  both  devices  in  a  single  circuit,  and  a  short  develop- 
ment of  this  appears  in  hopes  of  beginning  the  solution  to  this  prob- 
lem. 

The  analysis  is  to  be  done  with  the  thought  in  mind  of  obtaining 
a  set  of  design  equations  for  a  circuit  which  might  be  used  as  a  bi- 
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stable  network.  Both  the  two-branch.,  series  resonant  circuit,  which 
consists  of  series  elements  in  a  parallel  network,  and  the  two-mesh, 
parallel  resonant  circuit,  which  consists  of  parallel  elements  in  a 
series  network,  are  analyzed.  Consideration  of  triggering  schemes, 
although  of  great  importance,  occupies  a  small  portion  of  this  paper 
since  it  was  not  the  primary  purpose  to  design  an  actual  circuit  but 
rather  obtain  the  criteria  which  might  lead  to  a  realizable  circuit. 
The  production  of  an  actual  operating  circuit  awaits  the  development 
of  materials  with  lower  losses  and  slightly  greater  non-linearities. 
The  necessary  changes  and  the  progress  being  made  on  this  phase  of  the 
elements  are  entirely  beyond  the  scope  of  this  paper. 

The  author  wishes  to  express  his  deep  gratitude  to  Charles  A. 
Rosen  and  Charles  F.  Spitzer  of  the  General  Electric  Laboratories  for 
their  extensive  suggestions  and  guidance  in  the  completion  of  this 
paper.  Without  the  many  "skull"  sessions  with  these  gentlemen,  this 
paper  probably  would  not  have  been  written. 
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CHAPTER  I 
INTRODUCTION 

The  circuits  to  be  analyzed  in  this  paper  are  the  two-branch, 
series  resonant  circuit  shown  in  figure  I -la;  and  the  two-mesh,  para- 
llel resonant  circuit  shown  in  figure  I -lb.  Both  of  these  circuits 
utilize  only  one  non-linear  element.  In  chapter  IV  a  brief  develop- 
ment appears  of  the  two-branch,  series  resonant  circuit  using  two 
non-linear  elements  as  shown  in  figure  1-2. 

Inasmuch  as  these  circuits  contain  non-linear  elements,  methods 
other  than  linear  network  analysis  must  be  employed.  There  are  several 
approaches  for  analyzing  such  circuits  based  upon  developments  by 
Duffing,  Poincare  ,  Kryloff  and  Bogoliuboff  and  others  /T,  2,  3,  h, 
and  8J.     The  choice  of  one  approach  over  another  depends  on  the  accur- 
acy desired,  the  accuracy  obtainable,  and  the  exact  problem  confront- 
ing the  analyst.  The  method  of  Dryloff  and  Bogoliuboff  is  certainly 
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more  general  than  that  of  Poincarej  and  both  are  more  accurate  than 

that  of  Duffing.  On  the  other  hand,  the  method  of  Duffing  is  far 
easier  than  the  others  and  yet  yields  the  desired  information  to  a 
sufficient  degree  of  accuracy.  All  of  these  methods  mentioned  are 
analytical  and  suffer  from  one  inherent  disadvantage:  they  produce 
solutions  for  only  the  case  under  consideration  without  any  thought 
as  to  the  totality  of  solutions  available.  A  qualitative  method, 
known  as  the  topological  approach,  gives  an  overall  picture  of  the 
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Figure  I-la  Two  Branch,  Loss-Free  Series  Resonant  Circuit 
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Figure  I -lb  Two-Mesh,  Loss-Free  Parallel  Resonant  Circuit 


c, 


6 


Ls   sin  cut 
o 


41^U 


l 


■t.  = 


^ 


e»- 


-<:  - 


ii 
it 


Figure  1-2     A  two-branch  Loss -Free  Series  Resonant 
Circuit  using  Two  Non-linearities 


totality  of  solutions  but  yields  no  specific  information  on  any  partic- 
ular case  /~1>  3}   and  8_/. 

The  non-linear  capacitors  available  at  this  writing  had  relatively 
large  losses ,  were  dependent  on  temperature,  and  had  small  non-lineari- 
ties. Furthermore,  there  was  no  available  test  equipment  for  measur- 
ing the  component  parameters  at  the  higher  frequency  of  circuit  opera- 
tion. Finally,  the  determination  of  the  component  parameters  was  based 
on  some  rather  sweeping  assumptions  and  hence  accuracy  was  sacrificed 
ab  initio. 

For  these  reasons,  and  due  to  its  simplicity,  a  modified  form  of 
Buffing's  method  was  chosen.  In  order  to  clearly  understand  the  devel- 
opments of  chapters  II  through  IV,  a  brief  presentation  of  this  method 
and  the  reasoning  leading  to  it  is  given  in  what  follows. 

Consider  the  motion  of  a  simple  pendulum  having  no  losses  and  not 
under  the  influence  of  an  external  force.  The  differential  equation 
of  this  system,  as  shown  in  figure  1-3  is 

X    *  JL  <u^n   *  *   0  •  1-1 

The  solution  is  begun  by  using  the  substitutions 

Equation  1-1  then  becomes,  after  separating  the  variables: 


irJv  -    --f-A**>z  <Jx. 


1-2 


Integrating  this  expression  gives: 

2  J  1-3 

where  h  is  the  energy  constant  of  the  system.  From  this  equation  a 

plot  of  v  versus  x  is  made  as  shown  in  figure  I-lj.  From  equation  1-3 
and  the  plots  of  figure  I-lj.,  the  following  observations  can  be  made: 

1)  h  ">0  or  V*   becomes  negative,  which  defines  an  unreal- 
istic case: 

2)  if  ~Z2 £}\4.  £&  ,  the  curves  are  closed,  encircle  the  points 
v  =  0,  x  -  2nTV  ,  n  being  an  integer,  and  hence  the  pendulum  oscillates 
about  the  angle  x  -  2n1"f"  ;   and 

3)  if  h>^a  ,  v  is  never  zero,  and  thus  the  pendulum  con- 
tinues to  rotate  about  its  pivot. 

The  heavy  plot  corresponds  to  the  case  n  *  -£    ,  and  it  can  be 
shown  that  this  case  cannot  exist.  However,  this  value  of  h  divides 
the  totality  of  solutions:  those  that  are  periodic  and  those  that  are 
aperiodic . 

This  procedure  just  outlined  represents  the  crux  of  the  topologi- 
cal approach.  From  it  there  has  been  determined,  once  and  for  all, 
the  totality  of  solutions  in  the  system,  but  nothing  has  been  learned 
as  to  the  values  associated  with  any  particular  case. 

It  seems  evident  that  if  a  driving  force  were  included,  then 
whether  or  not  the  solutions  were  periodic  or  aperiodic  would  depend 
on  the  amount  of  energy  this  force  adds  to  the  original  system. 

Obviously,  if  aperiodic  solutions  are  allowed  in  the  electrical 
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Figure  1-3       Simple  Pendulum 
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Figure     I-I4. 


Velocity"  versus  displacement  curves 
for  simple  pendulum. 


circuit  cases,  then,  at  some  time,  the  currents  and/or  charges  will  be 
of  sufficient  magnitude  to  destroy  the  circuit  elements.  Thus  analysis 
can  immediately  refrain  from  these  solutions  and  concentrate  on  those 
in  which  periodicity  is  mandatory. 

In  order  to  obtain  these  periodic  solutions  we  now  consider 
Duffing's  method  which  was  first  applied  to  the  equation  bearing  his 
name: 

X+CZ    +&,*    H-Aa***    Fc#Lout.         x_h 

Note  that  the  function  of  x  appearing  in  this  equation  resembles  the 
first  two  terms  of  the  expansion  of  the  sin  x  function.  Hence,  it 
might  be  expected  that  the  solutions  are  much  the  same  as  those  of 
equation  1-1.  Such  is  found  to  be  the  case  within  the  region  of 
periodic  solutions  /~1  and  8_7« 

First,  assume  no  damping,  i.e.  c  -  0,  and  write  the  equation  in 
the  following  form: 

j>  =  -a,z  -aE  -xJ-+Fcooou±.        i-$ 

The  iteration  now  starts  with  the  first  approximation  (the  approxima- 
tion x  -  0  being  considered  trivial),  Xflsn^WWt"j  and  inserts  this 
into  the  right  hand  side  of  equation  1-5 •  This  obtains: 

Zl^-(^lk^0.iA3'F)t^U)t-  ±azk  Cta3out     1-6 
Integrating  this  expression  twice,  there  results 

the  constants  of  integration  having  vanished  in  order  to  preserve 


periodicity". 

At  this  point  Duffing  made  a  rather  far  reaching  choice.  In  lieu 
of  the  usual  procedure  of  prescribing  wand  then  determining  the  ampli- 
tude, A-,,  of  the  fundamental,  he  prescribed  A]_  -  A  and  thus  obtained 

LOl.r  0.,  +  f<L3A*-    J-    •  i_8 

It  is  noted  that  CU  becomes  a  single-valued  function  of  A  but  that 

A  could  be  a  triple-valued  function  of  UJ  .  It  was  for  this  reason 

that  A  was  prescribed,  and  hence,  UJ   determined. 

The  next  step  of  the  iteration  procedure  would  be  to  substitute 

the  solution 

X,*  Attest   4  ^A^c^3ujt 

into  the  right  hand  side  of  equation  1-5  and  again  integrate  twice 
letting  the  constants  of  integration  vanish.  This,  then,  is  Duffing's 
method.  For  rapid  convergence  a~  and  A  should  be  small. 

The  modified  procedure  used  in  this  paper  removes  the  restric- 
tions on  a_  and  A,  but  implies  that  the  driving  force  is  small  of 
order  a_  i.e.  we  set  F  =  a3^o« 

This  begins  by  adding  {jj*?(   to  both  sides  of  1-5  and  obtaining: 

%  +ufa  =  (w*-<UX  -a3"*3+  as6  tauft.  1-9 

Again,  iteration  starts  with  the  approximation  ^r  ActCCof  and  substitutes 
this  into  the  right  hand  side  of  equation  1-9. 
This  gives: 


1-10 
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Figure  I-5>a     A, 


Relationship  for  Doffing' a 
Equation 
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Figure  I -5  b     A 


Response  Curves  for 
Duffing ' s  Equation 


Now  to  perform  the  next  step  in  the  iteration  it  is  noted  that  the 
coefficient  of  Ctcuit=d  or  else  secular  terms  of  the  type    B"t  C6&  uft 
arise. 

This  procedure  gives: 

which  is  precisely  the  same  relation  resulting  before  if  a»F  is  re- 
placed by  F.  After  setting  the&Qujt  term  to  zero,  there  results: 

*,  +  u;1*,  =  -  -jr  (L&  A3  C^a  3  cot 

which  gives  as  its  solution 

*,=:  /4>(L*3CAjt  +  Jka^3  C#J3U>£.       1-12 
As  in  Buffing's  method  the  choice  is  to  make  Aj_  =  A.  The  procedure  is 
then  to  substitute  1-12  into  the  right  hand  side  of  equation  1-9  and 
repeat  the  steps  above.  The  iteration  is  continued  until  the  desired 
degree  of  accuracy  is  obtained. 

It  is  instructive,  and  will  be  useful  later,  to  plot  equation 
1-8.  Although  lo*  is  the  dependent  variable  in  this  case,  it  is  plot- 
ted on  the  abcissa  in  order  to  show  the  resemblance  of  the  response 
curves  to  the  linear  case.  The  plots  in  figure  I-£  are  for  the  loss- 
less case,  while  those  in  figure  1-6  are  for  the  case  where  losses  are 
considered  /  8  f*     It  might  be  added  parenthetically  at  this  point 
that  the  plots  in  these  figures  bend  to  the  right  or  left  according 
to  whether  a,,  is  greater  or  less  than  zero  respecitvely.  If  a~  is 
zero,  the  linear  case  results  and  the  curves  are  upright. 
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Figure  1-6 
Case 


Figure  1-7 
Capacitor 


Figure     1-8 


Critical  Point 


CjU 


for  Duffing !s  Equation  Low -Loss 


Q-E  Relationship  for  A  Non-Linear 


Nest  of  Q-E  Curves  for  Varying 
Field  Amplitudes 
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If  a  single  series  circuit  were  to  be  driven  with  a  constant 
frequency  source  whose  amplitude  could  be  made  to  vary,  then  operation 
would  take  place  along  a  vertical  line  such  as  that  shown  in  figure  I-£a. 
It  should  be  noted  that  as  the  amplitude  of  the  driving  force  is  slowly 
increased  a  point  is  reached  where  there  are  two  possible  values  of  the 
response,  A.  Since  A  represents  charge  in  the  electrical  case,  and 
current  is  the  time  rate  of  change  of  charge,  when  this  critical 
force  is  reached,  the  current  is  said  to  jump.  If,  instead  of  a 
single  series  circuit,  two  such  circuits  are  placed  in  parallel  as  in 
figure  I-la,  and  a  constant  current  source  used,  it  would  seem  plaus- 
ible that  when  the  critical  force  level  is  reached,  one  branch  current 
would  "jump"  up  while  the  other  branch  current  would  "jump"  down.  Such 
an  event  does  occur. 

The  question  now  arises  as  to  what  form  f(q),  shown  in  figure  1-1, 
shall  take.  If  a  charge  versus  voltage  plot  is  taken  on  one  of  these 
non-linear  capacitors,  a  curve  such  as  that  shown  in  figure  1-7  results. 
This  particular  curve  is  for  one  value  of  driving  voltage.  In  figure 
1-8  a  set  of  such  loops  is  shown  for  various  driving  amplitudes.  The 
heavy  line  connects  the  loop  tips.  Obviously,  no  simple  relation  de- 
fines the  Q-E  curves  of  either  figure  1-7  or  I-  8  - 

Since  the  real  interest  is  in  the  value  of  charge  which  can  cause 
bistability,  it  would  seem  evident  that  the  large  signal  loops  are  of 
more  importance.  This  is  due  to  the  fact  that  the  charges  must  be 
large  enough  to  cause  an  excursion  into  the  non-linear  region  at  the 
time  of  "jumping"  into  bistability. 
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If  a  single-valued  function  is  desired,  it  is  assumed  that  the  Q-E 
relationship  follows  the  mid-points  rather  than  the  upper  or  lower  por- 
tions of  the  loops.  This  tends  to  average  the  characteristics  and  has 
the  advantage  that  voltage  is  zero  when  the  charge  is  zero. 

The  assumptions  made  this  far  are  the  primary  limitation  on  the 
accuracy  of  any  solution.  It  would     seem  rather  irrational  to 
choose  an  involved  functional  relationship  however  accurate  it  may  be. 
Thus  the  defining  relation  is  chosen  as: 

e=  o.,^  +  <Lg  ^3i  1-13 

an  odd  function  which  is  required  for  the  midpoint  loci. 

The  necessary  information  has  now  been  obtained  in  order  to 
commence  the  analysis  of  the  bistability  problem. 
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CHAPTER  II 
TWO-BRANCH,   SERIES  RESONANT  CIRCUIT 

It  is  the  plan  of  this  chapter  to  investigate  the  two-branch, 
series  resonant  circuit.  Investigation  first  concerns  itself  with 
the  loss-free  circuit  in  order  to  determine  the  basic  relationships 
for  bistability.  After  these  conditions  have  been  studied,  emphasis 
will  shift  to  the  low-loss  circuit,  in  an  effort  to  determine  the 
effect  of  resistance  on  the  bi3tability  relationships. 


SECTION  I:       THE  LOSS-EREE  CIRCUIT 


The  equations  of  the  circuit  of  figure  I-la  can  be  written  as 
follows : 


Il-la 


e*  =  L  ££  +  ■  «,  +  bj&  ■ 

Subtracting  equation  Il-lb  from  Il-la,  there  is  obtained  one 
of  the  basic  equations  for  the  circuit: 

l  <JgJ£&>  +   e,-e,  -o.  u-2 

it 


The  other  basic  equation  for  the  circuit  is: 


?*-  l+p  ■ 


As  in  chapter  I  the  defining  relations  for  Ci  and  Q^   are 
chosen  as 


6,  =    a,  j,    +  a3  ^ 


H-lta 

Since  it  was  decided  that  the  circuit  should  have  a  constant 
current  source,  the  impedance  seen  by  the  generator  voltage,  S5  , 
must  be  that  offered  by  C^,  otherwise  the  varying  impedance  offered  by 
the  two  non-linear  branches  would  vitiate  the  constant  current  hypoth- 
esis.  Thus  is  leads  es  and  hence,  since  it  was  tacitly  assumed  that 
there  were  no  losses,  the  first  iteration  begins  with: 

Q,Q     =      A  XU^L-    Out  II-5a 

&0    a-      B  Jd^^V    6<jT  j  II-5b 

where  the  second  subscript  denotes  the  iteration  by  which  the  value 
was  obtained. 

Equation  II- 2  is  rewritten  as  follows: 

The  next  step  is  to  substitute  equations  II-5a  and  II-5>b  into 
the  right  hand  side  of  equation  II-6.  This  results  in  (see  appendix 
III  for  expansions  of  powers  of  trigonometric  functions): 
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ft'  >' +  ^l 'fa  -"  &>*-  %  )ft-8)- ^jt  g3| 


<kx*cot 


The  conditions  for  periodicity  require  that: 

Disregarding  the  trivial  case  of  A  =  B  and  using  equations 
II-3  and  II-5  there  results  two  equations  relating  A  and  B: 


%=0=    -S    + /m&ViAB+B1). 


II-10 


where 


S= 


UK  _  | 

a, 

j.  a-? 


ii-n 


11-12 


11-13 


&/    is  a  straight  line  whereas  ££  is  an  ellipse.  Using  the  rota- 


/ 


tion  matrix 

•V 

~£ 

-£1 

- 

B 

* 

A 

n-iU 


equation  11-10  becomes: 


11-15' 
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The  half  axes  are 

&«*  -  Y  ££n    )  and  11-16 

For  the  two-branch  series  resonant  circuit  the  second  iteration 
must  be  performed  in  order  to  obtain  bistability  conditions.  The  only 
information  obtainable  from  the  above  relations,  is  that  in  order  to 
have  real  solutions  of  y\2.  \    ^  ^*  0    and  hence   Uj  y>  GUo  where 
/  means:  "must  be  greater  than".  This  condition  can  be  deduced 
intuitively  as  follows .  For  bistability  A  •£?   B  and  hence  one  of  the 
two  identically  designed  branches  must  offer  less  opposition  to  charge 
movement  than  the  other,  due  to  the  non-linearity  of  its  branch  capaci- 
tance. This  requires  one  branch  to  be  "tuned"  nearer  to  CU  than  the 
other.  Now  for  a  capacitor  defined  by  an  equation  such  as  II -li,  it 
is  found  that  the  effective  capacitance  presented  to  the  fundamental 
voltage  decreases  as  the  amplitude  of  the  impressed,  fundamental 
voltage  increases.  Hence  as  the  amplitude  of  voltage  increases,  the 
resonant  frequency  of  the  branch  must  increase.  Hence  the  circuits 
must  originally  be  tuned  below  the  source  frequency  i.e.  CxJ>c<J0. 
Thus,  what  was  intuitively  believed  has  been  proven  by  the  inequality, 

&  $  o. 

Solving  equation  11-7^  after  the  condition  of  periodicity  has 
been  invoked,  there  is  thus  obtained 
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The  choice  as  explained  in  chapter  I,  is  to  make  A_  =A  and  B]_=B. 
Equating  the  portions  of  II -7  identifiable  with  q-,  and  q^  respectively, 
there  results 

0n    =    A  &>W   uft    -  hh      *Ur»   3WT!:  II-19a 


* 


,  -     Q/UsvyWt    —  DB    *Usn3oot  II-l?b 


where 


D  —      3      —     AVI 

"  SJLIO'L         ~      J*>(SH)--  n"2° 

Substituting  these  relationships  into  the  right  hand  side  of 
equation  II-6  there  results  (see  Appendix  II  for  details  of  the 
mathematical  development ) : 
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Invoking  the  condition  of  periodicity  and  disregarding  the  trivial 
on,  A  =  B,  thi 


solution,  A  =  B,  there  results: 


Using  the  rotation  matrix  of  equation  II-II4,  equation  11-22 
becomes 

+-£(&<,+3**'fsi'-hJio<y+0i).        11-23 

This  form  is  still  not  too  instructive.  However,  another  trans- 
formation reduces  the  form  to  where  the  desired  results  can  be  obtained. 
Letting  <?<-/>Ce<6      and  A   -  A4x^l&      ,  11-23  becomes: 

h  ■  ^  - .  *  +  * /fe^1  *  ^A  /^/      n-2^ 

where        J  fa)  r   J2  *  0d4  36 j  H-25a 

For  any  value  of  A  ^   has  only  two  finite,  real  roots  (posi- 
tive and  negative  values  of  the  same  magnitude)  and  hence  represents 
a  single  closed  loop.y^  (/^   and  (//   are  shown  in  figure  II -1.  The 
developments  of  these  curves  and  the  proof  of  the  single  root  magni- 
tudes are  shown  in  Appendix  II. 

From  figure  II-l,  the  mode  of  operation  may  be  explained.  For 
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small  values  of  Eg,  and  hence  Qs  ,  the  two  series  branches  behave  as 
linear  networks  and  thus  A  =  B.  As  Es  is  increased,  fy      moves  out  fr 
the  origin  parallel  to  the^  axis .  The  operating  point  remains  on  the 
intersection  of  $   and  the  line,  A  =  B.  In  Appendix  II  it  is  shown 
that  for  0-0  degrees,  the  slope  of  the  tangent  of  ^  is  parallel 
to  the  3   axis .  Hence  jfy   will  be  tangent  to  y£   at  the  operating  point 
Pj_.  When  yf  becomes  tangent  to  ^  ,  there  are  now,  two  other  operating 
points  available,  Pp  and  P  .  Although  no  mathematical  proof  is  given 
here,  it  can  be  deduced  from  the  work  of  Stoker  /~8"7  that  the  area 
around  Pj_  is  unstable.  Operation,  therefore,  "jumps"  to  either  P?  or 
P.,.  The  reasons  leading  to  this  deduction  of  P-^  being  unstable  are 
given  in  chapter  III  for  the  first  iteration  of  the  two-mesh  parallel 
circuit.  The  reason  for  postponing  this  discussion  until  then  is 
apparent  from  the  discussion  itself. 

A  circuit  in  which  the  controlled  switching  of  the  operating 
point  from  P  to  P~,or  vica  versa  is  obtainable  constitutes  a  trigger 

Cm  J 

pair. 

The  reader  might  now  infer  that  the  first  iteration  will  not 
predict  bistability  for  the  two-branch  series  resonant  circuit.  How- 
ever, in  chapter  III  it  will  be  shown  that  this  is  an  unwarranted  con- 
clusion; the  same  value  of  Es  which  will  cause  bistability  will  be 
predicted  by  two  different  approaches. 

For  operating  point  P-j_  in  figure  II-l,  the  following  values, 
developed  in  appendix  II,  are  obtained  (0  -  0°): 
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A    *  -V^JkL^S^^Jt^-aflTl^  n-26 

06  -  A+B***A    =  i$P     c  II_28 

-Ts  =   OJ  #s  n-29 

Obviously,  if  ty]   moves  out  from  the  origin  so  far  that  is  fails 
to  intersect  ^i    ,  bistability  fails  to  exist  and  both  branches  revert 
to  equal  currents.  Form  figure  II-l  this  appears  to  occur  when   is 
tangent  to  IfL   at  about  0  =  20°.  For  this  condition,  the  value  of  Ig 
becomes 

Greater  accuracy  can  be  obtained  by  finding  the  value  at  which 
vanishes,  but  this  hardly  seems  like  a  reasonable  procedure  in 
view  of  the  original  assumptions. 

It  should  be  added  at  this  point  that  after  bistability  occurs, 
there  is  no  assurance  that  increasing  Es  will  increase  Ig.  This  is  due 
to  the  phase  relationships  possible  between  A  and  B.  If  the  two  opera- 
ting points  are  located  in  the  first  quadrant  of  figure  II-l,  then  A 
and  B  are  in  phase.  If  the  two  operating  points  are  located  in  the 
second  or  fourth  quadrants,  then  A  and  B  are  180  degrees  out  of  phase. 
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SECTION  2:     THE  LOW-LOSS  CIRCUIT 


Attention  is  now  given  to  the  case  where  the  two  branches  contain 
losses  as  shown  in  figure  II-2.  Equation  I 1-6  is  rewritten  to  include 
the  effect  of  resistance  as  follows: 

Assuming  that  R  is  small  compared  to  Xc-, }   then  Is  will  still  be 
primarily  determined  by  the  impedance  of  C-j_.  For  the  remainder  of  the 
discussion  in  this  chapter,  it  will  be  assumed  that  k)>y  B.  Hence,  I, 
is  very  nearly  in  phase  with  Is.  Since  Is  will  lead  Es  by  90  degrees, 
the  first  iteration  begins  with  the  approximations: 

§/c    -      A    <U^.   uot  II-32a 

fa*   B<^(*>*  -9*)  IJ"32b 

Substituting  these  values  into  the  right  hand  side  of  equation 
11-31,  there  is  obtained  thereby: 


-[(^-%)P-5fl3BiJ<W(^-^) 
_L  .2.3  43  A^Uji.  ~   &3  B*A^3(ujt-0), 

*  hl  «l  rJ  n-3: 
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Figure  II-2 
Two-Branch  Low-Loss  Series  Resonant  Circuit 
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Figure  II -3 

Phase  Relationships  for  Two-Branch,  Low-Loss  Series 
Res  onant  Circuit  * 

2k 


In  order  to  preserve  periodicity  all  of  the  fundamental  component 
terms  of  LU  must  vanish  identically.  Hence  there  results: 


a. 


Expanding  the  terms  involving  (U/t-p)   and  equating  the  coefficients 
of  sin  wt  and  cos  wt  respectively,  there  results: 

[-SB+^BB]C*5p  +  ^BA^y*    L-SA   4/mA3J   II-35a 

The  phase  diagram  appears  in  figure  II-3*  Since  R  was  presumed 
small,  ^  is  very  nearly  zero  degrees  and  hence  cos  J> £?  1.  Using 
this  fact  and  dividing  equation  II-35>a  by  II-35b,  there  results: 


or 


_s+^(A>^B^Bi)=-^fi.s^r^  n-36 

If  R  is  allowed  to  vanish,  the  first  iteration,  loss-  free  case  is 
immediately  derivable  from  the  above.  Experimentally  Ar\  g  is  of  the 
order  8  x  10"-^  whereas  s  is  unity  or  greater,  so  that  the  assumption^ 
S>>mS)"is  not  too  unreasonable.  If  it  is  allowed  that  s  »tvy\B*  then 
11-36  reduces  to 
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'a  /rvi    u»(5/Wi         11-37 

Now  ^  will  not  be  a  real  illipse,  and  hence  bistability  cannot 


72 

result,  unless 


R  <•    ^> 


CO  '  11-38 

which  places  an  upper  limit  on  the  losses. 
For  &  %  0        11-37  gives 


11-39 


If  R  is  allowed  to  vanish,  the  first  iteration,  loss-free  case 
is  obtained,  as  can  be  seen  from  equation  11-10. 
From  figure  H-2,  it  can  be  seen  that: 


or 


£o  =  ^.[tsAt/Vh^A^ut  +  qpAcHvtl.   ii- 


10- 

The  magnitude  of  E0  can  be  found  by  first  finding  its  square  as: 

Substituting  the  value  of  A  from  11-39 ,  and  upon  simplification, 
there  results: 


F»-"*K£-#&Y/+^).    1T-h3 
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Since  6  <^  Q 


?**  ?**  + 1  *>* -&y 


=  y®?r+^ 


and,   since 

T      a    U)  B    -  £° 

R-     kiSVV—  -   ^^^     ) 

p-    ^Fl^i     ^s/w  y  , 


n-w 


n-l*5 


The  ratio  of  A  to  B  (i.e.  the  ratio  of  the  larger  to  the  smaller 
charge  magnitudes)  is  defined  as  the  switching  ratio, A   ,  and  can  be 
obtained  from  equations  11-39  and  11-1*5: 

Now  from  equations  II-lja  and  Il-ljb,  remembering  that  B^O,  there 
results: 

f(   =   &,A  f"  ^^  II-i,7a 

£  x  =   °n  B  .  II-ltfb 

Since  the  usefulness  of  this  circuit  is  predicated  upon  a  suffic- 
iently large  difference  between  the  two  branch  output  voltages,  it  may- 
be desirable  to  solve  for  this  voltage  difference:   Hence,  by  sub- 
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traction,  and  by  using  equations  11-39  and  II-2i5  there  results: 

which  for  the  loss-free  case  gives 

Thus  it  can  be  seen  that  although  the  switching  ratio, A   ,  is 
independent  of  the  degree  of  non-linearity,  m,  that  the  voltage  diff- 
erence is  inversely  proportional  to  its  square  root.  It  should  be 
noted,  however,  that  A   is  not  the  ratio  of  E-,  to  E?,  for  the  circuit 
under  consider ati on «  This  value,  still  independent  of  m,  is  given  by: 


=  A  [i  +  s-  £). 


11-50 


The  smaller  the  value  of  m  is,  the  larger  the  potential  difference 
between  the  non-linear  capacitors  is  for  a  given  ^  .  This  follows 
from  equation  II-1|9  and  the  fact  that  A  is  related  to  the  parameter, 
s,  by  equation  II-I4.6.  This  potential  difference,  however,  is  obtain- 
able only  by  making  A  larger  than  B.  nut  since  the  switching  ratio, 
d\   ,  has  been  prescribed,  both  A  and  B  increase,  with  B  remaining  in 
the  linear  region  and  A  in  the  non-linear  region  of  the  Q-E  relation- 
ship. 

However,  although  A    is  independent  of  m,  it  must  be  remembered 
that  too  small  a  non-linearity  results  in  currents  of  such  magnitude 
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that  they  may  destroy  the  components  before  testability  is  reached. 

All  of  the  above  relations  were  derived  from  the  first  iteration 
of  the  circuit  with  small  losses.  The  second  iteration  is  now  com- 
menced in  order  to  find  what  modifications,  if  any,  need  to  be  made 
to  the  first  iteration  relations. 

After  the  periodicity  conditions  have  been  invoked,  equation 
11-33  becomes 

whose  solution  is 

),  -/^  -  (A,  -B,)*u^Cot  -  DA7*^  3*ot 

-hbB<2*r?3(io£  -j^>)  ,   H-52 
Again  choosing  A-j_=A  and  B]_=B  there  results,  for  the  second  approxima- 
tions : 

0„   «="  /)**nu>£-  bA^un  Scot         n^3a 

Substituting  these  expressions  into  the  right  hand  side  of  11-31, 
gives: 


f 


*•   V 
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-t 


«■        »  ♦       /     /    I      ' 


where  the  dots  refer  to  terms  involving  ,?£*»£,  SuSC^k/t  <£r>cJ S^Ujt*t 

The  conditions  for  periodicity  become: 

-  ^  (A-B)  +  ^  D B*^,  3y  =  O  ii-5S 

-  |«3  b^-B*)]  -     O  .  XI-56 


From  11-55 ,   the  phase  relationship  is  determined  as 


/  3  a3        &> 

$>  i 

Thus,   since  <u*n~£>f   —     '  . 

K   -   TuJ    (A-B  ' 


II-S8 


Solving  for  A-B  results  in: 


Assuming  that  equation  II-lj.6  does  in  fact  give  a  reasonable 
approximation  forz^  ,  there  results,  upon  substituting  into  11-59 

.  s 

A   —    S-/mg^  '         n-60 
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Here  again  since  bistability  implies  a  value  of  ^  greater  than 
unity,  a  small  value  of  m  may  necessitate  a  value  of  B  so  large  that 
the  components  themselves  might  be  destroyed. 

Substituting  equation  11-57  in  equating  11-56  and  by  using  the 
first  two  terms  of  the  expansion  of  Ceo.3  0 3   there  results,  assuming 
(/?   to  be  small: 

v      '  11-61 

Collecting  terms,  excluding  the  trivial  solution,  A=B;  applying  the 
rotating  matrix,  II-lljj  and  using  the  polar  coordinate  form,  there 
results  for  0  =  0°: 

Except  for  large  values  of  R  and  ^P  ,  l/~  is  approximately  equal 
to  AJL  as  given  by  equation  11-27. 

Briefly  summarizing  these  results: 

1)  in  order  to  have  bistability 

Is  rr  -&£>(.  l5t^A6244*JM--#J%   *       H-26 

where     ^  =r  "#^W "7  &'  '& S  + *8°5      \  H-27 

2)  the  switching  ratio  is  given  by 

^  uTR       j 
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3)  circuit  losses  should  be  such  that 

R  < 

h)   the  difference  of  non-linear  capacitor  potentials  is  given  by 


CO    )  H-39 


E-E^^o+^yx 


^v)  J  11-1*9 

5)  the  ratio  of  the  two  non-linear  capacitor  voltages  is: 

|^=    a(i+S  -  i»);  W         n-so 

6)  the  non-linearity  should  be  such  that 

e   S  g,  ->    A    >  I  n-60 

s-  />^S^ 

does  not  require  B  so  large  that  the  component  ratings  are  exceeded. 
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CHAPTER  III 
TWO-MESH,  PARALLEL  RESONANT  CIRCUIT 

The  plan  of  this  chapter  is  to  develop  relations  for  the  two- 
mesh,  parallel  resonant  circuit  similar  to  those  discussed  in  chapter 
II  for  the  two-branch  series  resonant  circuit.  It  will  also  be  shown 
that  the  requirement  of  a  constant-voltage  source  in  the  place  of  the 
previously  used  constant-current  source  gives  rise  to  a  different  solu- 
tion curve  relating,  A,B,  and  Es,  the  source  voltage. 

As  in  chapter  II,  the  first  consideration  will  be  the  loss-free 
circuit  and  then  the  low-loss  circuit. 


SECTION  1:   THE  LOSS-FREE  CIRCUIT 


The  equations  for  the  circuit  of  figure  I-lb  can  be  written  as 
follows: 


J,   =      l£a£  +    4&  m-ia 


I    _      j  re*  &  +    sfe.  ■  i"-lb 


Subtracting  equation  Ill-lb  from  III -la  there  is  thus  obtained 
one  of  the  basic  equations  for  the  circuit: 
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-tfa-etft   +       <£$>'&   -  0.  UI-2a 

Since  there  are  no  d.c.  sources  in  the  circuit,  no  loss  of 
generality  occurs  if  equation  III-2a  is  differentiated  with  respect 
to  time.  This  differentiation  results  in 

This  equation  is  the  same  as  was  obtained  in  chapter  II  for  the 
two  branch,  series  resonant  circuit. 

The  other  basic  equation  of  the  circuit  is: 

<?/  +  ea   -  es  .  in-3 

As  with  the  series  circuit,  the  defining  relations  for  <3,  and  CA 
are: 

<?,   =    °->  ^  +  4sf?>  ni-U 

In  lieu  of  the  constant  current  source,  a  constant  voltage  source 
is  desired.  Hence  0,  should  be  fairly  large.  Thus,  i0  is  a  small 
quantity.  For  lossless  circuits,  &,     and  ^  both  lag  iQ  by  90  degrees 
if  U)  >  CUo  .  The  charge  will  be  in  phase  with  the  voltage  and  the 


first  approximations  can  be  written  as: 

l0    =      A  COO   COt)  ln_5a 

-     B   CflO   U)t.  III-£b 


% 


Proceeding  as  in  chapter  II,  the  two  solution  curves  relating  A 
and  B  are: 
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fesO=  -  S  +  /m(Ai'-hAB+  Bx)  ,  in-6 

It  should  be  noted  that  the  $£  curve  is  common  to  both  systems,  but 
that  the  straight  line,  ff  ,  has  now  become  the  cubic  ^  .  These 
curves  are  shown  in  figure  III-l,  as  derived  in  appendix  II.  The  fact 
that  %   is  a  cubic  as  shown,  makes  it  possible  to  demonstrate  bistabil- 
ity  from  the  first  iteration.  The  operating  point  for  low  values  of 
Es  is  on  the  intersection  of  fa   and  the  line,  AsB.  As  E  is  increased 
(ft,  moves  out  from  the/^  axis,  its  asymptotic  line  remaining  parallel 
to  theS  axis.  When  ^becomes  tangent  to  ^  at  P]_,  the  operating  point 
"jumps"  from  P,  to  P?  or  P_.  Mathematical  proof  that  P.  is  an  unstable 
operating  point  is  beyond  the  scope  of  this  thesis.  As  stated  in 
chapter  II,  however,  the  instability  may  be  inferred  from  the  work  of 
Stoker  [_  8_/.    The  reasons  leading  to  this  conclusion  of  instability 
will  now  be  shown. 

In  chapter  II  it  was  shown  that  the  half  axes  of  the  ellipse  are 


<L~Y\ 


3  s»n      j     and 


III-8 


3*  =  /^pT-  1Z1'9 

The  first  of  these  represents  the  "distance"  from  the  origin  to 
P-i  .  Since  A=B  at  this  point: 

Substituting  this  value  into  equation  III-7  results  in 
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es=  ao.,o+  4)yi^; .  m-ii 

As  suggested  in  chapter  II  there  is  another  method  by  which 
equation  111-10  and  hence  equation  III-ll  may  be  derived.  It  is  the 
criterion  of  this  method  which  leads  to  the  conclusion  that,  P,  is 
unstable.  In  chapter  I  it  was  shown  that  the  frequency  for  the  loss- 
free  Duffing  equation  was: 


uj2  =  a,  -4-  ^a%A 


'■3#     ~  ~A      '  III -12 


Since  Duffing's  equation  is  the  same  as  the  equation  for  one  mesh 
of  the  circuit  under  consideration  if  F  is  replaced  by  the  voltage,  E; 
across  the  mesh,  the  frequency  is  given  by: 

w»,  A,   +  §.*sAl-  JL.  in-13 

^    ¥  i  A 

The  locus  of  the  points  where  the  curve  of  Fig.  I-3>b  has  vertical 
tangents  can  be  obtained  by  differentiation  of  equation  III -13  with 
respect  to  A  and  setting  QUO /a A  equal  to  zero.     This  procedure  results 
in: 

602-%    -|«3   A*  =    0  m-ii, 

The  value  of  A,  for  a  particular  LU ,   which  lies  on  the  locus 
given  by  III-lli  is  found  by  solving  equations  111-13  and  III-llj.  simul- 
taneously. This  results  in  the  value  of  A, 


n  ~  r  3/m     *  iii-io 

Stoker  has  shown  that  the  region  bounded  by  equation  111-lij.  is 
unstable  /  8  7°  It  was  assumed  that  both  meshes  of  the  two-mesh  cir- 
cuit were  identically  designed,  and  that  the  circuit  was  supplied  from 
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a  constant  voltage  source.  Thus,  when  A  reaches  the  value  given  by 
111-10,  only  one  of  the  mesh  voltages  can  "jump"  up.  The  other  "jumps" 
down.  Hence,  P-^  is  unstable  and  operation  "jumps"  to  either  Pp  or  P^. 
The  value  of  Es  given  by  equation  III-ll  has  another  interpreta- 
tion which  may  be  obtained  by  study  of  figure  III-l.  If  E  were  larger 
than  the  value  given  by  equation  III-ll,  then  there  is  the  possibility 
of  having  two  sets  of  bistable  operating  points,  since  ^  will  inter- 
sect ^  in  four  points.  Thus  for  no  more  than  two  possible  stable  oper- 
ating points 

Es  <:aa,6+^)7££"-  111-15 

This  apparent  contradiction  between  equations  111-15  and  III-ll  is 
explained  as  follows.  Bistability  results  when  Es  reaches  the  value 
given  by  equation  III-ll.  However,  as  the  amplitude  of  Is  is  increased, 
there  is  no  assurance  that  the  value  of  Es  will  increase  due  to  the 
phase  relationships  possible  between  A  and  B.  As  stated  in  chapter 
II,  operating  points  lying  in  the  first  quadrant  represent  the  in  phase 
conditions,  whereas  those  operating  points  lying  in  the  second  or 
fourth  quadrants  represent  the  180  degree  phase  conditions.  However, if 
Es  can  be  caused  to  exceed  the  value  given  by  equation  111-12  after 
bistability  has  been  reached,  then  four  stable  points  may  be  possible. 

As  Es  is  further  increased,  there  is  a  value  beyond  which  bistabil- 
ity will  cease  and  reversion  to  the  A=B  condition  occurs.  This  value 
of  Es  is  determined  by  finding  the  common  point  of  tangency  between  £z 
and  #  besides  the  one  already  found  at  Pj_  .  For  ^_and  M  to  osculate 
at  a  given  point,  the  following  relations  must  be  satisfied  at  the  point 
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"Ml      V&, 

"\A  ^A  m'16 

><A.  _   >J&  III17 

when  Es  has  been  adjusted  so  that  the  point  of  tangency  lies  on^_, 
then  both  the  necessary  and  sufficient  conditions  have  been  met. 
Applying  equations  III-16  and  111-17  to  equations  III- 6  and 
III-7  there  results,  respectively: 

/+  ZsmA7'  =    smfaj+S),  ni-18 


/-f  3/*n£Z  =    "*"  (28 +A). 


gives 


III -19 
Mvision  of  these  two  equations  by  each  other,  and  simplification 

3/9n(A-B3)+  (6*»AB-/)(A-8)  =  o*  m-20 

The  solution  A=B  is  trivial  for  the  present  development  and  hence 


there  results: 

to=0   =  A^+ZAB+B*   -  3?m 


111-21 


Application  of  the  rotation  matrix 


A 
B 


'//3   -  '//a 


I/. 


to  equation  111-21  results  in 
a  hyperbola  with  half -axes. 


3/^n- 


ys>yri 


4*=^ 


3sm 


111-22 


III-23 

III-21* 

111-25 
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The  ellipse,  (^  and  the  hyperbola,  %-  will  intersect  if  the  half 
axes  obey  the  relations 

From  this  inequality  another  criterion  for  bi stability  results: 

S   >  5^  J         or  111-27 

Uf    >  jr  Ufi.  HI-28 

The  two  criteria  for  bistability,  i.e.  the  required  values  of  s 

and  Es,  thus  far  derived  have  been  obtained  from  the  first  iteration. 

The  second  iteration  should  give  corrections  to  these  relations. 

Proceeding  as  in  chapter  II  the  two  solution  curves  thus  obtained  are: 

^3  =  0  =    -£  +  (A^+AB  +3  z)  +tfA  U*J3+AZB  Z+A3*+  B*) 
+  3&Ytf+A%AlB+#  &+£&•*&+&)  ni-29 

The  curves  are  shown  in  figure  III-2  Continuing  the  reasoning  of 
the  first  iteration  procedure,  bistability  occurs  when^j^  becomes  tan- 
gent to  y^   at  Pi  .  Applying  the  rotation  matrix^  equation  111-22  to 
equation  111-29  and  equation  111-30  gives,  respectively: 


Uo 


ia 


As  before,   the  solution  curves  are  now  brought  into  the  polar 
form  through  the  equations :  <=*  y&C&o  &     ax\&£=jO&^\@    .      This  results  in: 

%S6=  -£  ^fty^tyU^^'  m.33 


where,  f(0),  g(0),  and  h(0)  are  as  in  chapter  II;  and 

/M  -  30*0.0   -   aC6c36  IH-35a 

k(6)  -  S"C**4    -  VC4*S~6  ni-35b 

/(&)   -  7l*4  6+'VC*636  -3tc*d*Z  +  gCe*>4.  lll-35c 

Figure  III-2  shows  a  plot  of  these  curves.  The  curve  y£   differs 

but  little  from  %~~}   except  for  large  values  of/*??  or /&  .  The  curve, 
U?   shown  is  therefore  an  approximation. 

The  value  of  Es  required  for  bistability  is  obtained  by  substitu- 
ting into  equation  111-31;  the  values:  @-  O  degrees  and 


111-36 


where 


jc    =    *^-'(*'<iJ*«r).  ni-37 


This  procedure  results  in 

E»  =  i  &>  *  "7?  *+  *&/9 r  +  ^r /?  7)>       to-38 

where  x>    is  given  by  equation  III -36, 

There  is  no  easily  derived  correction  for  the  required  value  of  s. 
Application  of  equations  111-16  and  III -17  to  equations  111-29  and 
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111-30  results  in  an  equation  which  is  very  tedious  to  solve.  Con- 
sidering the  approximations  first  made  regarding  the  Q-E  relationship 
of  the  non-linear  capacitor,  the  solution  does  not  seem  worthwhile. 
SECTION  2:  THE  LOW-LOSS  CIRCUIT 

As  in  chapter  II,  attention  is  now  given  to  the  low-loss  circuit 
as  shown  in  figure  III-3,  with  the  phase  relationships  shown  in  figure 
III-1|.  Equation  III-2b  then  is  rewritten  as: 

($-£)+«>%-£)•  **(frfr)-k(i-i)-tte-**)-  m-39 

Assuming  R  to  be  large  compared  to  Xc-, ,  the  first  iteration  commences 
with  the  approximation: 

$*=  Ao*>  €^t  in-ta 

ae  =r   B  t*a(to?+  S)  in-iiOb 

Substituting  these  into  the  right  hand  side  of  equation  111-39, 
and  involving  the  conditions  for  periodicity  there  results: 

Using  the  expansions  of  cos(^-f£)  and  sin  (tot-tS)   and  equating 
coefficients  of  like  functions  there  results,  since  cos^Ai.  /: 

Dividing  equation  III-l|2a  by  III-i|2b  and  subsequent  simplification 
results  in: 


1*3 


6 

Is  sin  cut 
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Figure  III-3:    Low-Loss,  Two  Mesh  Parallel  Resonant  Circuit 
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Figure  III-I4: 


Phase  Relationship  for  Low-loss,  Two-mesh 
Parallel  Resonant  Circuit. 
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Since  in  practice,  />w£  is  of  the  order  of  8  x  10" ^   it  can  be 
neglected  in  comparison  with  unity.  Thus  equation  III-li3  reduces  to 


2     2 

Solving  this  equation  for  A  -frAB-fB   results  in: 

A*+AB+BX=  s  ,  + 2F 


in-16 


Clearly  this  reduces  to  the  loss-free  case  as  R  becomes  infinite. 
If  the  right  hand  side  of  equation  III~i;5  becomes  negative,  the  equa- 
tion fails  to  define  a  real  ellipse.  Thus  in  order  to  permit  bistabil- 
ity,  the  following  inequality  must  hold: 

For  B %d       j  A  is  taken  from  Ill-itf?  as: 

Solving  equation  III -la  for  i05  and  using  the  fundamental  com- 
ponent only  results,  after  simplification,  in  the  expression: 
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Since  B  is  small,  the  admittance  of  the  associated  mesh  is: 


^         S3-  .     III-li9 


Also,  since  B  is  small,  there  results  from  equation  III-l|b: 
E^  =  0.  , B   =    Ta  ;   or 


III-50a 


Again,  defining  the  switching  ratio  .,Z\  ,  as  the  ratio  of  A  to  B 
i.e.  the  ratio  of  the  larger  to  the  smaller  charge  amplitudes,  there 
results: 


__£     (1^-ZEk) 


^  =  OuL    S  -h  J  111-51 

Now,  from  equations  III-lja  and  Ill-l+b,  and  setting  J3&0   ,   there 
results : 

£j   =•  O^B-  III-52b 

Since  the  usefulness  of  this  circuit  is  predicated  upon  a  suffici- 
ently large  difference  between  the  two  mesh  voltages,  it  may  be  desir- 
able to  solve  for  this  voltage  difference.  Hence,  by  subtraction  and 
use  of  equations  III-U7  and  III-£Ob,  there  results: 


(Mx      "I     r       .,/.„!  A.  <&X 


r     r   =  >Vy    ■    Ja,  [j-  Mfe^L  +  flj£LZUlll-53 
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For  the  loss-free  case,  this  reduces  to: 

Here  again,  although  the  switching  ratio  is  independent  of  the 
degree  of  non-linearity,  /yr>)   the  voltage  difference  is  inversely  pro- 
portional to  its  square  root. 

The  ratio  of  the  mesh  voltages  is  given  by 

I -*(■ +  *£$)■      — 

Since  the  second  iteration  of  the  low-loss  circuit  produces 
algebraic  equations  which  are  not  readily  solvable  and  the  first, 
iteration  of  the  low-loss  circuit  has  become  rather  involved,  it 
hardly  seems  reasonable  to  attempt  a  second  iteration  in  an  effort 
to  obtain  more  information. 

In  summary,  then,  the  following  conclusions  apply  to  the  two- 
mesh  parallel  circuit: 

1)  The  driving  frequency  must  be  such  that  for  bistability  to  be 
possible, 


cu5  >  a)  r  j  m-28 

2)  The  minimum  value  of  voltage  required  for  bistability  is, 
from  the  first  iteration, 


:>e  such  that 

QJL  .  I II -I16 


-s 
3)  The  losses  must  be  such  that 


R  > 


is   > 


hi 


k)   The  switching  ratio  is      , 

-   /  ^  _  uolL   \ 


) 


iii-£l 


$)   The  voltage  difference  is  given  by 


ill -53 


iteration 


where 


and 


6)  The  value  of  Es  required  for  bistability  is,  from  the  second 


^^^+Wfr+*&Y) 


f^<L  -*.9/7] 


A-    t^^(^^) 
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CHAPTER  IV 


USE  OF  TWO  NON-LINEAR  ELEMENTS: 
FIRST  ITERATION 


In  chapters  II  and  III,  the  two  basic  trigger-pair  circuits 
utilizing  a  non-linear  capacitor  have  been  briefly  studied.  The  two- 
mesh  series  resonant  circuit,  using  a  non-linear  inductance,  has  been 
analyzed  by  Spitzer  /  7__7"«  Intuitively  it  is  felt  that  the  inclusion 
of  both  types  of  non-linearities  would  reduce  the  currents  required 
for  the  same  switching  ratio.  This  chapter  will  investigate  this 
possibility  by  analyzing  the  loss-free,  two-branch  series  resonant 
circuit  using  two  non-linearities  as  shown  in  figure  1-2.  Due  to  the 
complexity  of  obtaining  a  solution  and  the  inherent  inaccuracies  pro- 
duced by  the  assumed  current  and  charge  functions,  the  analysis  will 
be  confined  to  the  first  iteration. 

For  reasons  given  in  chapter  I,  the  current- flux  and  charge- 
voltage  relationships,  are  chosen  as: 

/(<£)  SA  -  -^  5^  +4*  &  IV-1 


IV-2 


Furthermore,  since  ^ 

-d"  <«",  IV-3 


/ 


equation  IV-2  can  be  written  as: 
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IV-8b 


One  equation  for  the  system  is  given  by 

whereas  the  other  equation  is: 

<,       +  u£      =    ^  IV-6 

The  branch  currents  and  voltages  are  from  equations  IV-1  and 
IV-2  respectively:* 

U,  =  ^7  <&,    •+  ^3   <&/  iv-7a 

^a  =  ^  ^    /  -4  ^/  I*-7b 

65    =    ^    V-  0*  g?  IV-8a 

S   =    **£    *  *s/ 

Again,  the  capacitor,  C^,  is  to  be  the  primary  impedance  offered 
to  the  driving  voltage.  Hence  Is  leads  Es  by  °0  degrees.  For  the  loss- 
free  circuit,  the  branch  currents  will  be  in  phase  with  the  source 
current,  Is*  Thus  the  branch  fluxes  will  also  be  in  phase  with  Is, 

and  the  first  approximation  is  taken  as: 

c/>,4     =    /I  C#<2>  Uft  IV_0a 

<2^     =    S  C*2    6t/7?t  IV_ob 

Equation  IV-5  can  be  differentiated  with  respect  to  time  without 
affecting  the  solutions  since  there  are  no  d.c.   sources.     After  per- 
forming this  differentiation,   the  equation  is  solved  for  (<&,  —  <j£>2    ). 
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In  the  resulting  equation,  tv* (<&-<&$  is  added  to  both  sides,  and  there 
is  thus  obtained: 

The  iteration  begins  by  substituting  the  approximations,  <^/(9  and^ 
into  the  right  hand  side  of  equation  IV-10  and  then  invoking  the  condi- 
tions for  periodicity. 

The  most  expedient  way  to  do  this  is  to  write  <2.  as 

e  -  a,  djf  *  az  ^f  . 

The  value  of  <?  is  then  obtained  from  equation  IV-3  and  the  result 
cubed.  The  value  of  <?-*  is  then  differentiated  directly.  This  proce- 
dure eliminates  the  tedious  progress  of  computing  9  and  then  multi- 
plying  by  4  . 

From  equation  IV-3,  using  equations  IV- ?a  and  IV- 9b  there  results: 

+  &,  (A'-Bl)*~>  set  ■  IV"12 


Letting 


and 


F=    4L(A-8)+^(A*-B3)~]  lv.13a 

£  =   ^  "(^-B3);  IV-l3b 


there  results  upon  cubing  IV-12: 


at  ^  2(£L?  +  F(^)^u^cut    -l(F*<r+<r/0***tf' 
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,3 


—  -&    xl<c*n    P&jlt. 


IV-lli 


The  condition  for  periodicity  hence  becomes 

.3 


The  terms  F  and  G  are  replaced  by  equations  IV-13a  and  IV-13b 
respecitvely  and  upon  simplification  there  results: 

-% (A-B) -h  (/+4p)(A*-S?)+  ztm (A*- 8V 
+  &fJ«n«(A1-B->)+KJL„HA'*-BC)=d    IV.16 

where 

ATI  =r    ^  &  IV-17a 

^     —     %  Z?  IV-l?b 

J    =     -&-  IV-l?c 

^\ 
S    P*       7%a  -/  IV-17d 


2,    _i  a,i, 


IV-17f 


U)*  -      -_L-  _        ^-»^i  IV-17g 


The  result  of  substituting  equations  IV-9a  and  IV- 9b  into  IV-6  isi 

Is  =  J-,  (A  +  B  )  +  ^  ^  V*  5J.        iv.i8 
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This  is  the  equation  for  ^  found  in  chapter  III.  Omitting  the 
trivial  solution,  A  -  B,  equation  IV-16  defines  a  solution  curve,  tyi    . 
Bistability  results  when^  becomes  tangent  to  (/C    .  This  occurs  at  a 
point  on  the  A  =  B  line.  To  find  the  value  of  Is  required  for  bi- 
stability, equation  IV-16  is  first  transformed  into  the  polar  coordin- 
ate form  by  the  same  transformations  used  in  chapters  II  and  III.  The 
resulting  equation  is  solved  f  or/2  ,  and  this  value  is  used  to  compute 
A.  Equation  IV-18  then  produces  the  required  Is. 

The  transformations  indicated  above  reduce  equation  IV-16  to 

+&<A"«*'y+&M«^  ivi9 

where  l/^l3^)   and  ^(&)    are  defined  as  before,  and 

jft)m  &&4  *6  -tfbafe  +i?c*i  g»  +3&fo  H.       iv-20 

o 

If  the  term,  £> }     is  neglected,  equation  IV-1°  is  of  the  same  form 
as  the  equation  defining  pi   .  Since  the  functions  involving  the  co- 
ordinate $   ,  are  always  positive,  the  effect  of  the  second  non-linear- 
ity is  to  reduce  the  required  current  for  bistability.  The  amount  of 
reduction  depends  on  the  value  of  l,m,  and  n. 

There  is  no  explicit  statement  derived,  for  the  value  which  s  must 

take.  However,  from  previous  developments  of  the  single  non-linearity 

M 

circuit,  it  can  be  assumed  thatS>0>  at  least. 

Substituting  the  value  0  =  0  degrees  into  equation  IV-19  results 
in: 
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+  ffijm^t+jft  ^,„y*.  Qi        I¥.20 

This  equation  is  a  quartic  in  /O     .  The  general  form  of  the 
quartic  is  written  as 

X*+  clx3  +  4x*+  a^c  +</=  a,  IV_21 

where  a,b,c  and  d  are  real,  constant  coefficients. 
To  find  the  roots  of  equation  IV- 21,  it  is  necessary  to  find  any- 
real  root  of  the  cubic: 

%  3-*^  Wd£-*^  -[cVfff«'.  9i  J7  =  (7.       17-22 

The  roots  of  equation  IV-21  are  then  taken  from  the  two  quad- 
ratics: 

x  V  [f  -  Yf +*>,-<*]  *+(%,  +  >Jy^7>a    ""23a 
1 V  [S  +/f*^J*  *&  -i^T^c   IV"23b 

where  ^  is  the  root  found  in  equation  IV-22. 

Obviously  this  is  a  tedious  process.  Fortunately,  it  needs  to 
be  done  for  only  one  value  of  0.  However,  to  arrive  at  a  general 
solution  in  terms  of  the  parameters  l,m,n,  and  s  is  not  within  the 
realm  of  this  thesis.  For  any  given  components  and  source  frequency, 
the  value  of/O  obtained  through  the  quadric  solution  is  substituted 
into  the  polar  form  of  equation  IV-17: 
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This  procedure  gives  the  required  value  of  I_  for  bistability. 
In  summary,  the  results  of  analysis  of  the  doubly  non-linear 
series-resonant  trigger  pair  are  as  follows: 

1)  The  value  of  Is  should  be  less  for  bistability  in  the  double 
non-linearity  circuit  than  for  the  single  non-linearity  circuit; 

2)  There  is  no  explicit  development  for  the  value  which  s  should 
take,  but  it  is  implied  by  previous  study  that  S  ~>  0. 

3)  The  required  value  of  Is  is  obtained  from  auxiliary  quartic, 
cubic  and  quadratic  equations,  as  given  by  equations  IV-20,  IV-21  and 
IV-22. 
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CHAPTER  V 
EXPERIMENTAL  RESULTS 

The  experimental  results  obtained  with  the  circuits  under  consid- 
eration are  given  in  this  chapter.  They  are  not  offered  especially  to 
confirm  theoretical  developments,  but  to  demonstrate  that  the  circuits 
are  operable.  If  the  experiments  agree  with  the  theoretical  develop- 
ments, so  much  the  better. 

The  charge-voltage  relationships  were  determined  through  the  use 
of  the  test  equipment  arrangement  shown  in  figure  V-l  for  60  cycle 
supplyj  and  by  the  arrangement  shown  in  figure  V-2  for  the  5  to  20 
kilocycle  range.  Sample  hysteresigrams  for  frequencies  of  60,  5000, 
10,000  and  20,000  cycles  are  shown  in  figure  V-3.  In  figure  V-lj  is 
the  hysteresigram  taken  at  60  cycles,  but  with  the  vertical  gain  and 
intensity  set  so  that  the  loop  characteristics  are  distinguishable. 
The  locus  of  the  tips  of  a  nest  of  20,000  cycle  major  hysteresis 
loops  is  shown  in  figure  V-5. 

For  all  experimental  work  the  MUCCN  type  VSE  non-linear  capaci- 
tors were  used.  The  values  of  the  coefficients  for  these  type  capaci- 
tors were  obtained  from  figure  V-5  by  approximating  the  average  initial 
and  final  slopes.  The  Q-E  relationship  resulting  from  this  procedure 
is  shown  as  the  dashed  curve  in  figure  V-5« 

The  coefficient  values  thus  obtained  are 

Ol  ,    r   9*  M  volts/coul  v~la 
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Figure  V-l    Test  Equipment  for  obtaining  hysteresigram  at 
60  cycles.  The  VTVM  should  be  peak  reading. 


Figure  V-2   Test  equipment  for  obtaining  hysteresigrams  from 
5  kcs  to  20  kcs.  The  VTVM  should  be  peak  reading, 
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f :  20  kcs 
Emax:  350  volts 
Qmax:  0.27  ucoul. 


f :  10  kcs 
Emax:  350  volts 
Qmax:  0.3  ucoul, 


f :   5  kcs 

Emax:  300  volts 

Qmax:  0.299  ucoul. 


f :  60  cycles 
Emax:  375  volts 
Qmax:  0.315  ucoul. 
(Q  axis  reversed) 


Figure  V-3.   Hysteresises  of  MUC0N  type  VSE  non-linear  capacitors. 


Figure  V-4-.   Hysteresigram  of  MUC0N  type  VSE  non-linear  capacitor 
showing  open  loop  for  a  source  frequency  of  60  cycles. 
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*9 


a3  =   9./6a/o*' **tf/c4tJs 


V-lb 


V-lc 


The  circuit  of  figure  I-la  was  first  tested  using  the  following 
component  values: 


PEST  NO. 

1 

2 

3 

L 

110  rah 

1.62  mh 

0.1*6  rah 

F 

20  kcs 

233     kcs 

233     ks 

For  these  components,  the  computed  and  experimental  results  are 
tabulated: 


nssT 

1 

2 

3 

Out 

8.l8xl09 

S.^xlO11 

1.96xl012 

s 

0.9h 

2.86 

.091 

Is   (comp) 

73  ma 

1.59a 

k&h  ma 

Is   (exp) 

70  ma 

.&lia 

327  ma 

Since  the  parameters,  a.  and  a^,  were  not  separately  determined 
for  the  higher  frequency,  it  is  not  suprising  that  the  experimental 
and  theoretical  results  disagree.  However,  from  tests  two  and  three, 
it  was  noted  that  the  higher  value  of  s  gave  a  quicker  "jump"  but  that 
the  temperature  of  the  crystal  oven  in  which  the  capacitors  were 
located  rose  10  degrees  centigrade  suggesting  a  radical  change  of  the 
capacitor  parameters  as  a  function  of  temperature. 
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Figure  V-6.  Ej_  and  E2  for  the  two-branch  series  resonant  circuit. 
The  source  voltage  is  larger  for  the  waveforms  on  the  right. 


Figure  V-7.  E^  and  E2  for  the  two-mesh  parallel  resonant  circuit. 
The  source  current, Iif  is  larger  for  the  waveforms  on  the  right. 


WOllr        f 


Figure  V-8.  Relative  switching  ratios  and  switching  rates  of  the 
single  non-linear  resonant  oircuit  and  the  double  non-linear  cir- 
cuit. Figure  a  shows  the  switch  for  the  two-branch  aeries  resonant 
circuit  using  a  non-linear  inductance;  figure  b  shows  the  switch  for 
the  two-branch  eerie*  resonant  circuit  using  both  non-linear  induct- 
ance and  non-linear  capacitance. 
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11111111111 
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1000  pps 


300  pps 


100  pps 


Figure  V-9.  Output  voltsge  for  vsrious  trigger  rates  of  one  branch 
of  the  two-branch  series  resonant  circuit  utilizing  two  non-linear 
elements.  The  voltage  waveform  is  taken  across  one  of  the  non-linear 
capacitors.  The  trigger  pulse  was  6.1  microseconds  in  duration  and 
65  volts  in  amplitude. 
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Switching  was  not  obtainable  by  use  of  a  pulse  generator,  but 
only  by  the  use  of  a  long-time  (five  seconds  or  more)  short  circuit 
across  the  high  current  capacitor.  This  indicates  that  the  capacitors 
were  heating  excessively  and  that  the  subsequent  change  in  parameter 
values  was  detrimental  to  switching.  Methods  of  maintaining  both 
capacitors  at  the  same  constant  temperature  were  not  feasible  because 
of  limitation  of  time  available.  Figure  V-6  shows  the  output  voltages 
E,  and  Ep  at  233  kcs  after  the  bistability  condition  has  been  obtained. 
The  phase  relation  conditions  should  be  noted.  The  in-phase  condition 
was  first  obtained,  and  by  increasing  Es,  the  out-of -phase  condition 
resulted. 

The  parallel  circuit  was  next  investigated  at  the  frequency  of 
20  kilocycles.  The  data  on  this  test  are  given  below. 


TEST  NO. 

k 

L 

110  mh 

*>? 

9.1&X109 

S 

0.91* 

Es   (comp) 

1*82  volts 

Es   (exp) 

578  volts 

Figure  V-7  shows  the  output  voltages,  E,  and  E2  after  the  bi- 
stability condition  has  been  obtained.  Note  that  the  two  voltages 
begin  out-of -phase  and  approach  the  in-phase  condition  as  Is  is 
increased. 

Losses  in  both  the  series  resonant  circuit  and  parallel  resonant 
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circuit  cause  the  large  deviation  from  the  180  degree  relationship. 

The  fifth  test  concerned  itself  with  .he  double  non-linear,  two- 
branch  series  resonant  circuit  versus  a  single  non-linear  two-branch 
series  resonant  circuit.  The  noii-linear  inductors  were  those  used  by 
Spitzer  /~7_7»  The  parameters  of  these  were: 

/uy  =■  2  x  10^  amps/oeo 
/  *    2   x  1017  amps/aO)3 


/n  =   7-5  x  1012  c&Cx 


V-2a 
V-2b 

V-2c 
V-2d 


N  =   200 

The  capacitor  replacing  the  non-linear  capacitance  for  the 
single  non-linearity  circuit  was  chosen  to  give  approximately  the 
same  switching  ratio, 

The  results  are  given  below: 


Test  $ 

Nonlinear  Capacitor 
Linear  Inductor 

Nonlinear  Capacitor 
Nonlinear  Inductor 

F 

233  kcs 

233  kcs 

I« 

h92   ma 

hlS  ma 

Figures  V-8a  and  V-8b  show  the  relative  switching  ratios  and 
switching  rates  of  the  two  circuits.  The  double  non-linear  circuit  is 
shown  to  switch  in  about  two  cycles  of  the  supply  voltage  whereas  the 
circuit  using  only  the  non-linear  inductance  switches  in  about  two  and 
one-half  cycles. 

The  sixth  test  consisted  of  operation  of  the  double  non-linear 
two-branch  series  resonant  circuit  at  various  switching  rates.  The 
results  of  this  are  shown  in  figure  V-°.  A  Hewlitt-Packard  Model  212A 
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Pulse  generator  was  used  as  the  trigger  source. 

Only  during  a  brief  interval  of  time  v  as  switching  obtained  on 
the  two-branch  series  circuit  using  a  non-linear  capacitor.  It  was 
apparent  from  the  experiments  that  the  temperature  change  of  the  capac- 
itors is  the  primary  deterrent  to  switching. 

It  is  felt  that  the  data  presented  confirms  the  theory  of  the 
first  four  chapters  to  a  good  degree.  The  design  equations  seem  to  be 
valid  for  operation  near  the  frequency  region  in  which  the  capacitor 
parameters  were  determined.  At  high  frequencies,  the  losses  are  so 
great  in  comparison  to  the  degree  of  non-linearity,  that  the  circuit 
does  not  operate  within  the  theory  propounded  by  the  analysis. 
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APPENDIX  I 
TRIGGERING  SCHEMES 


It  was  stated  in  chapters  II  and  III  that  a  circuit,  in  which 
controlled  switching  is  obtainable  between  the  two  bistable  points, 
would  constitute  a  trigger  pair.  Obviously  the  circuit  must  also 
include  means  for  obtaining  this  triggering.  The  purpose  of  this 
appendix  is  to  present  a  few  ideas  as  to  how  this  may  be  obtained. 
The  discussion  will  be  limited  to  the  suggested  methods  without  too 
much  thought  as  to  their  output  requirements.  One  important  question 
is  left  unanswered  viz.  by  what  method  is  memory  produced.  This  is 
an  important  question,  but  is  considered  beyond  the  scope  of  this  paper. 

When  the  pendulum  was  studied  in  chapter  I,  it  was  found  that  the 
division  between  periodic  and  aperiodic  solutions  was  dependent  on  the 
energy  of  the  system.  The  circuits  studied  in  chapters  I,  II,  and  III 
were  shown  to  pass  from  the  unistable  condition  to  the  bistable  condi- 
tion  when  a  particular  value  of  Es  or  Is  was  obtained.  Although  the 
implication  was  not  forthright  the  approximation  of  the  Q-E  functions 
of  these  circuits  to  the  sine  function  used  in  the  pendulum  case  leads 
to  the  conclusion  that  the  two  states -unistable  and  bistable -are  con- 
cerned with  the  amount  of  energy  possessed  by  the  individual  branches 
or  meshes. 

As  an  approximation,  only  the  energy  contributed  by  the  funda- 
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mental  component  will  be  considered. 

Consider  first,  the  two-branch  series  resonant  circuit.  There 
are  two  possible  methods  of  trigger  application:  a  pulse  applied  to 
the  capacitor j  or  a  pulse  induced  in  the  coil.  In  figure  AI-1,  both 
methods  are  shown. 

From  equation  II-375  the  value  of  A  for  which  bistability  occurs 
is  given  by 


A  =  fz 


/TTi  OL^S/m  AI-1 

and  from  equation  11-39  the  value  of  A  for  the  maximum  A    in  the  bi- 
stable condition  is  given  by: 

rim*}     y  <*n        et^s/m  M~2 

The  current  in  one  branch  is  ^^  and  for  each  of  the  two  states- 
unistable  and  bistable  -  it  is  given  by: 

T,      =   ^A  M-3a 

Since  the  energy  is  associated  with  the  magnetic  field  for  part 
of  the  cycle  of  the  driving  frequency,  the  energy  contained  by  the 
branch  can  be  written: 


[J  *  $  LJ/  —    /-ft  Aj^a 


for  the  unistable  state, or 


LU* =    ^  L  A""*>t  ^b 


for  the  bistable  state. 
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Figure  AI-1 


Two  Trigger  Schemes  For  Two  Branch 
Series  Resonant  Circuit. 
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Figure  AI-2 


Two  Trigger  Schemes  For  Two-Mesh 
Parallel  Resonant  Circuit 
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The  energy  difference  is  thus  given  by: 

Um#*    U  -      3/rn   V  CL?  S     J  AI-$a 

by  using  equations  AI-1  and  AI-2.  This  is  the  amount  of  energy  that 
this  branch  must  lose  before  it  is  possible  to  switch  it  into  the  lower 
charge  state. 

Since  it  was  assumed  that  B  <&6  ,   the  amount  of  energy  this  cir- 
cuit must  gain  by  virtue  of  a  pulse  is  given  from  equation  AI-lia  as 

When  branch  "1"  has  lost  the  amount  of  energy  specified  by  equa- 
tion AI-3>a  and  branch  M2W  has  simultaneously  gained  the  amount  speci- 
fied by  equation  AI-6a,  switching  can  occur.  In  order  to  assure 
switching,  the  equality  should  be  replaced  with  the  ^inequality  sign 
as  follows 


AI-5b 


Since  the  trigger  pulse  is  applied  to  both  circuits,  it  must 
cause  this  exchange  of  energy  to  be  as  specified.  At  this  point  the  . 
investigation  should  inquire  about  memory.  Although  the  subject 
memory  is  beyond  the  scope  of  this  paper,  a  short  attempt  will  be  made 
to  explain  how  the  energy  exchange  may  be  caused. 

The  trigger  pulse  adds  energy  to  both  branches,  "l"  and  "2".  The 
addition  of  energy  to  branch  "l"  does  not  change  the  effective  capaci- 
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tance  of  the  non-linear  element  by  much  since  it  was  already  operating 
in  the  saturated  region.  However,  the  energy  addition  to  branch  H2M 
does  cause  a  large  change  in  the  effecting  capacitance  of  the  non-lin- 
ear element.  Hence  the  opposition  to  the  fundamental  component  de- 
creases to  the  point  of  non-linear  resonance  in  branch  "2"  but  does 
not  decrease  appreciably  in  branch  "l" .  Since  the  source  is  provid- 
ing a  constant  current,  the  current  now  begins  to  divide  between 
branches  °1"  and  "2",  adding  to  branch  "2"  by  taking  from  branch  nl". 
Branch  "2"  successively  offers  less  opposition  to  the  current  whereas 
branch  "1"  successively  increases  its  opposition.  The  process  cascades 
until  switching  has  occurred.  When  the  energy  in  branch  "1"  has 
reached  the  value  specified  by  equation  Al-lja,  the  trigger  pulse  should 
cease.  Otherwise  it  would  be  possible  to  have  branch  "1"  return  to  its 
original  state  of  heavy  conduction.  The  reason  for  this  is  that  the 
pulse  is  supplying  a  certain  fixed  amount  of  energy  which  must  divide 
between  the  two  branches.  Since  the  effective  capacity  of  the  unsat- 
urated branch  is  greater  than  that  of  the  saturated  branch,  the  unsat- 
urated branch  is  more  receptive  to  the  energy.  However,  when  the 
branches  reach  the  energy  level  specified  by  equating  AI-iia,  the  recep- 
tion is  equally  well  for  both  branches .  Only  the  effect  of  cascading 
can  then  cause  the  switch. 

The  least  amount  of  energy  that  the  pulse  must  give  is  that  speci- 
fied by  equation  AI-6a.  Since  the  energy  will  divide,  and  at  some  time, 
equally,  perhaps  the  pulse  energy  required  should  be: 
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The  required  energy  expression  for  the  two-mesh  parallel  reso- 
nant circuit  can  be  similarly  obtained. 
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APPENDIX  II 


DERIVATION  OF  THE  SECOND  ITERATION  SOLUTION  CURVES 


The  purpose  of  this  appendix  is  to  show  how  the  second  iteration 
solution  curves  were  obtained.  The  second  iteration  begins  with  the 
approximations : 

0lt  =    A**™^  'fid  *^  5ur£  AII-la 

02l    =  B }  <Us>>iC*)t-DB*U<ri3C(jt  .  All-lb 

These  are  substituted  into  the  right  hand  side  of  equation  II-6,  which 
after  omission  of  the  trivial  solution,  A=B,  results  in: 

&=0   =  -S+/7»foA3+6*)*tn&$fS4ffrA&4!) 


+Sm/)Y^3W*£hf3s+MAB&J.k11-' 


First,  the  linear  transformation: 

A 
B 


AII-3 


is  applied.  This  results  in 

This  equation  is  next  transformed  into  polar  coordinate  form 
by  the  relations  cSj£>  &**&  0     and  j0-=^>&t*?  0    .  This  results  in 
(see  appendix  III  for  expansion  formulas): 
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where 


^)  =  X5>  J  &«■?£  -  0,  s-Ota  *y  AII.6b 

/( #,)  =  /*"  *  b>cMy-  /ec*>2&  -ce&6}+.  7c„fo     AII.6c 

Equation  AII-5  is  now  reduced  to  a  cubic  through  the  relation 
ps/l>     This  results  in: 

From  equations  AII-6a,  AII-6b,  and  AH-6c,  it  is  noted  that, 
for  all  6}  ^(&)  >  OjQ(&)>0  and  J(j)  >  O  .  Thus  the  terms  of 
equation  AII-7  varying  in  sign  only  as  does  J\* 

The  coefficient  signs  are 

%(*):     —     +      +       + 

0C-4):    -    -    +     ~  • 

Thus  there  is  one  or  no  possible  positive  roots j  and  there  are 
two  or  no  possible  negative  roots.  Since  there  are  three  roots,  it 
can  be  stated  that  there  are  two  negative  roots  and  one  positive  root. 
The  interest  lies  entirely  with  the  positive  root  since  the  negative 
roots  produce  complex  values  of  /2   which  have  no  meaning  as  far  as  the 
circuit  is  concerned.  The  "value"  /)=©£>  does  not  satisfy  the  equation. 
Therefore  ^  is  a  single  closed  curve. 

It  now  becomes  necessary  to  find  the  real  root  of  r '-^/» 
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The  general  form  of  a  cubic  equation  is 

&xs  +sA  y.1-^  sc  *  */=o.  AII_0 


The  trigonometric  method  of  solving  cubics  is  used  and  the 
real  root  is  given  by 


i  """'"" 


g£ ;  AII-9a 


where 


'^/'^^f   6*6  L*  ^    'FfJ  AII-9b 


or 


U    «i  <5     and        f+*3  ~    °J 
i.Zf<0        and      tfV  Xx>0-j 


A II-  9c 


AII-9d 


and        ^  -    #£  —  ^  *  AII-9e 

/  =  4(30.e4-CL*c/)-s4*  AII-9f 

The  upper  sign  is  chosen  in  the  J/,  evaluation  if  A>0   and  the 
lower  sign  if  n*0.     It  turns  out  that  for  all  &j  9>0   and  thus 
equation  A II- 9a  is  the  vehicle  by  which  the  roots  of  ^  are  obtained. 

The  equations  of  &(fL.)     are  tabulated  below  for  values  of  0  . 

0  %   (r  ) 

0°  -s/m     |  1.5  p     41  -2^D  r2   41-75  D2r3 

10°  -s/m    4  l.ltfr     +1.0  D  r2   -fl.86  D2r3 

20°  -s/m     41. 385 r    J1.23D  r2   42.065D2r3 

30°  -s/m     41.25   r    41.19  Dr2   42.22  D2r3 
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0 

^   (r) 

1*0° 

-s/m 

+1.09  r  h 

^  1.07    D  r2+  2.15  D2r3 

If 

-s/m 

+       r       \ 

f.  1.00    D  r24     2     D2  r^ 

50° 

-s/ra 

-f0.913r  • 

[•  0.905  D  r2+  2.07  D2**3 

60° 

-s/m 

40.75  r  - 

f  0.687  D  r2-f     1.28D2r3 

70° 

-s/m 

|0.6l7r  - 

\  0.+70  D  r2+  0.753D2r3 

80° 

-s/m 

fO.53  r  i 

F  0.31D  r2     +  0.37  D2r3 

90° 

-s/m 

+0.5   r    ^ 

f  0.251)  r2     +  0.25  D  r 

From  these  equations,  the  value  ofyO  is  obtained  for  each  angle 
specified 
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K      .88ZS+.687 


US*  +  .tot r  *  Tx      •****  +  ' 
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I 


.^  ....  v* 


\ 


kC 


\  *  I 


These  are  the  only  values  needed  to  form  a  plot,  since  ^3  is 
symmetrical  about  both  the  <=*£.  axis  and/^  axis  as  can  be  seen  from 
equation  AII-1*. 

The  plots  in  figures  II -1  and  III-2  were  obtained  from  the 
above  tabulations  for  an  s  value  of  0.68.  Although  not  labeled, 
the  coordinates  are  in  units  of  m  2. 

The  development  of  %   and  y£   is  done  in  a  similar  fashion. 
Tabulated  equations  of  m  ,  for  values  of  0  are  given  below.  Since  $£ 
differs  from  j£  but  little,  except  for  large  values  of  X>   or  m,  it 
has  not  been  tabulated.  The  only  real  interest  lies  in  the  fact 
that  both  fy   and  W  have  only  one  real  positive  root.  For  the  value 
0  =  90  degrees,  all  coefficients  of  Z)  vanish.  Hence,  /y   and/^  are 
open  curves.  Since  they  were  derived  from  monotonically  decreasing 
functions,  ^i(9)  and  £a.(<?)  by  a  linear  substitution,  they  will 
also  be  monotonically  decreasing  functions  about  the  axis  of  symmetry. 
Hence  the  shape  is  as  shown  in  figures  III-l  and  III-2. 

6  % 
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3d* 


^) 


~jZ£si  +  At£3f>    +/S/y&3 


Using  the  method  previously  described  for  finding  roots  of  the 
cubic,  there  arise  solutions  such  as 


*,  =  i,WfiW"/(fi4l, 


A 11-10 
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where  b  and  c  are  constants  determined  for  each  value  of  0.  For 
values  of      n —   ^  I  >   the  approximation 

^0^~^=Jn  U+  y^V/  )  Z&  ^n  (Z+/)  AII-11 

can  be  used.  Since  the  approximation  specifies  X^)>  there  is  the 
series  equality*. 

^*y/)«  x-ixV^xV......       Ani2 

Using  only  the  first  term  of  AII-12,  equation  AII-10  becomes 

Again  since    -  yT   £<\  ,   the  hyperbolic  sine  function  can 
be  replaced  by  the  first  term  of  its  power  series  expansion.  Thus 

*->   ^  36L,  All-lfc 

For  the  values  of  m  and  a-,  used  in  experiments,  Es  should  be 
less  than  100  volts  for  all  the  approximations  to  be  valid.  It  was 
on  this  basis  that  7y  and  yjk   were  plotted  in  figures  III-l  and 
III-2  respectively.  It  is  not  anticipated  that  a  marked  difference 
would  result  if  the  value  of  Es  were  such  that  equation  AII-10  must 
be  used. 

The  slope  of  the  fy  curve  is  obtained  by  application  of  the 
formula 
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where  T  is  the  angle  the  tangent  makes  with  the  horizontal  axis. 
Since  the  horizontal  axis  is  the  oL  axis,  for  the  tangent  to  be 
parallel  to  the  .^  axis,  the  denominator  of  AII-15>  must  vanish,  or 
the  numerator  must  become  infinite. 

Differentiating  A II- 5  with  respect  to  9,  there  results: 

From  equations  AII-6a,  AII-6b  and  AII-6c,  it  can  be  seen  that 
f'(9)j  g'(9)  and  h'(0)  each  have  the  factor  sin  0,  sin  29,  or  sin 
40.  Thus,  /O    vanishes  for  9-0  degrees  and  therefore  so  does  the 
denominator  of  AII-15*  It  has  then  been  proven  that  the  tangent  of 
^L   at  P_  is  parallel  to  the   axis. 
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APPENDIX  III 
MATHEMATICAL  FORMULAS  USED  IN  ANALYSIS 


This  appendix  lists,  for  the  readers'  reference,  mathematical 
formulas  used  in  the  analytical  portion  of  this  paper. 

The  trigonometric  expansions  are  obtained  by  applying  the  bi- 
nominal expansion  to  the  exponential  form  of  the  sine  and  cosine  func- 
tions. Thus  there  results 

(a^t/yix)     =       \       2  ■  UII-la 

Alll-lb 


AIII-2a 


AIII-2b 
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^y\/T)X) 


£srn+l 


(-')  \.AC+*\frn\+h*fi)mX+>n 


AIII-3 


(Cm,  njf^L     ^L  f(L*<Lfa#) 


/n  Z 


AIII-li 


In  determining  that  y^was  an  ellipse  and  investigating  if  it 
became  an  imaginary  ellipse  due  to  the  presence  of  losses,  the  follow- 
ing properties  of  conic  sections  were  used. 

The  general  form  of  a  conic  is  written  as 

The  relations,  D,  d  and  c  are  defined  as: 


D  = 


a 

h 

g 

h 

b 

f 

g 

f 

c 

AIII-6a 


d  = 


a  h 
h  b 


AIII-6b 
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e  =  a  +  b  AIII-6c 

The  conic  is: 

1.  a  parabola  if  d  :  0,  B\  Oj 

2.  parallel  lines  (may  be  coincident  or  imaginary) if 
d»0,  D=0; 

3.  an  ellipse  if  d>0,  eD<0; 

U.  no  locus  (imaginary  ellipse)  if  d  >  0,  eD>0; 
5>.  a  point  ellipse  if  d>0,  D—0; 
6.  an  hyperbola  if  d<  0,  D^  0 
7«  two  intersecting  lines  if  d<"0,  D="  0. 
In  order  to  remove  the  xy  term  of  the  generalized  conic  equation, 
the  rotation  matrix 


L/«J 


is  applied,  where  0  is  given  by 
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